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The Ginzburg-Landau (GL) equations for a d x i_ y i superconductor are derived within the context 
of two microscopic lattice models used to describe the cuprates: the extended Hubbard model and 
the Antiferromagnetic-van Hove model. Both models have pairing on nearest-neighbour links, con- 
sistent with theories for d-wave superconductivity mediated by spin fluctuations. Analytical results 
obtained for the extended Hubbard model at low electron densities and weak-coupling are compared 
to results reported previously for a d-wave superconductor in the continuum. The variation of the 
coefficients in the GL equations with carrier density, temperature, and coupling constants are cal- 
culated numerically for both models. The relative importance of anisotropic higher-order terms in 
the GL free energy is investigated, and the implications for experimental observations of the vortex 
Oh! lattice are considered. 
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I. INTRODUCTION 

There is mounting experimental evidence to sjiggest that the high-temperature cuprate superconductors have an 
' order parameter with unconventional symmetryoB Indeed, recentplosephson interference measurementsQ are strongly 
indicative of an order parameter with d x 2_ y 2 (d-wave) symmetryH which has line nodes along \k x \ = \k y \. The linear 
density of states associated with the resulting low-energy excitations is thoughtO to account for the linear temperature 
\ dependence of the specific heatO as well as the linear temperatureEl and magnetic fieldu dependence of the penetration 
depth found for YBa 2 Cu 3 7 _a (123land Bi 2 Sr 2 CaCu 2 O a (2212). 

A number of experimental resultsa are consistent only with an order parameter of combined s-wave and d-wave 
symmetry. Sigrist and RiceEHl have shown that in weakly orthorhombic cuprates (such as 2212 or 123) a small s-wave 
component would be present in addition to a critical d-wave order parameter. In tetragonal systems (such as the thal- 
lium compounds) thai-favour d-wave superconductivity, however, an s-pssjawe component can only bp-.nuclcated locally 
near inhomogeneitieaid such as domain walls, impurities|L3 or vortices. E3ii3 Indeed, Soininen et al.i3 investigated the 
structure of an isolated vortex for a d-wave superconductor within Bogoliubov-De Gennes theory and found that a 
non-zero s-wave component is induced in the vortex coca. Their results, interpreted within the context of the relevant 
phenomenological Ginzburg-Landau (GL) free energyJIH imply a non-trivial topological structure for the additional 
s-wave component .E-3 As a consequence, the supercurrent and magnetic field distributions for an isolated vortex near 
H c i exhibit a fourfold anisotropy in proportion to the magnitude of the s-wave component.c3 In addition, the vortex- 
lattice structure near Zf c2 deviates significantly from the usual triangular Abrikosov lattice, becoming increasingly 
oblique with increasing s-wave admixture I 14 " 16 ! 

It remains uncertain, however, whether the anisotropy in the structures of an isolated vortex and the vortex lattice 
is indeed predominantly due to the admixture of an s-wave component. Ichioka et al&3 have shown that the inclusion 
of higher-order d-wave gradients in the GL free energy can give rise to a fourfold-symmetric current distribution 
around a vortex even in the absence of an induced s-wave component. One of the objectives of the present work is 
the clarification of this issue. Indeed, we find that the contributions to anisotropy of a fourth-order gradient term and 
the s-wave component are comparable, and tend to compete. 

While phenomenological GL theory has been highly successful in predicting many interesting properties of d-wave 
superconductors in external fields, the relative magnitudes of the various coefficients appearing iu-thc free energy and 
their dependence on temperature, filling, and field are presently unknown. An earlier derivationEB of the free energy 
from a continuum model could not include lattice effects that are believed to be important in theories of d-wave 
superconductivity. Consequently, as will be discussed later, certain technical difficulties arose which would not have 
appeared in the continuum limit of an appropriate lattice model. In any case, it would be useful to derive the GL free 
energy using models relevant to the high-T c oxides. In the present work, the GL equations are derived microscopically 
within the context of two such models: the extended Hubbard and 'Antiferromagnetic-van Hove' models. 

The extended Hubbard (EH) model, which includes a nearest-neighbour attraction in addition to the usual on- 
site repulsion, is one of the simplest lattice models which allows for a d-wave superconducting instability. Pairing 
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occurs along nearest-neighbour links, appropriate for theories where d-wave superconductivity is mediated py an- 
tiferromagHfttic-fluctuations (see Ref. || and references therein). It has been employed in several arpJyticalE] and 
numericall3'E2lEil investigations of d-wave superconductivity. The EH model has recently been shown,E3 however, to 
favour d-wave superconductivity only in a very small parameter space, preferring a phase separated or spin-density 
wave state. i_. 

The 'Antiferromagnetic-van Hove' (AvH) mode£3 strongly-jjavours d-wave super.GQnductivity while incorporating the 
coexisting antiferromagnetic correlations observed in NMR,cJ neutron scattering,E3 and angle-resolved photoemission 
spectroscopy (ARPES)Ej experiments. High transition temperatures are obtained in the model due to the presence 
of a van Hove singularity in the hole density of states near the Fermi energy. An extended, flat band near (§,■§) in 
momentum .soace is consistent with numerical investigations of a single hole propagating through an antifcrromagnetic 
backgroundEZI, and with experimental evidence.E3 In the AvH model, holes are constrained to move within a single 
sublattice of a uniform antiferromagnetic background in order to minimize frustration. The hopping parameters are 
chosen to best fit the quasiparticle dispersion for YBCO measured using ARPESra 

In Section II, the Ginzburg-Landau equations for the gap functions and supercurrent are derived microscopically 
for both the EH and AvH lattice models using a finite-temperature Green function method. The relations defining 
the transition temperatures are investigated in Section III. It is found that only a d-wave transition is favoured for 
the AvH model; the d-wave transition temperature ~ 100 K is consistent with the high-temperature (high-T c ) 
oxides. The EH model, in contrast, can have either an s-wave or d-wave instability. S-wave is favoured at low electron 
densities while d-wave is favoured either at high densities or at lower densities with strong on-site repulsion. The 
equations for T s (the s-wave transition temperature) and Td are found analytically in the limit of weak-coupling and 
low electron densities. The corresponding analytical solutions for the AvH model are difficult to obtain due to the 
complicated angular dependence of the AvH dispersion. The GL free energy is derived for both models in Section IV. 
The coefficients of the GL equations are found analytically for the EH model in the same limit described above. The 
coefficients are calculated numerically for the EH model and for the AvH model near optimal doping. In section V, 
we summarize our results and discuss the experimental implications of the GL equations we have derived. 

II. THE LATTICE GL EQUATIONS 

The Hamiltonians for the extended Hubbard (EH) and antiferromagnetic-van Hove (AvH) models are respectively: 

(ij)a it? i (ij)aa' 



H^= tll £ dc^fi^+tzo £ ctc^^-^rw-^ J2 n ia n jr ,„ (2.2) 

{(ij))a («y)»<7 ia (ij)<?<?' 

where = c[ a c ial A is the vector potential associated with the external magnetic field, <j) = hc/e is the flux quantum, 
and fj, is the chemical potential included to fix the density. In the EH model the carriers are electrons, and positive Vq 
and Vij imply on-site repulsion and nearest-neighbour attraction, respectively. The superconducting carriers in the 
AvH model are holes propagating through the antifcrromagnetic background of the undoped 'parent' state. Second 
and third nearest-neighbour hopping parameters are respectively t\\ — 0.04125 eV and £20 = 0.02175 eV.a The 
absence of nearest-neighbour hopping in the AvH model reflects the restricted Hilbert space of the carriers; holes are 
constrained to move within a single spin sublattice in order to minimize frustration and preserve antiferromagnetic 
correlations. The values of t\i and £20 are chosen to result in a large density of states near the bottom of the hole 
band, located at (ir/2,ir/2) in momentum space. The coefficient of the nearest-neighbour attraction V = 0.075 eV 
is chosen to yield a d-wave transition temperature Td ~ 100 K at optimal doping (/1 rj —0.075 eV, or hole density 
(n) ~ 0.2), consistent with the high-T c oxides.E^I 
If the lattice sites i and j are nearest-neighbours, we can write the mean-field EH Hamiltonian 

A* H 4(r)c CT (r) + IK ( r ) c l ( r ) c T ( r ) + H - c ] 

r,<r r 

-A* s (r)c r (r) Cl (r + $)+H.c], (2.3) 



H™(B) =-t^ 4(r + S)c a (r)e^ - 

r,8,a 

-I^[A|(r) Ci (r) CT (r + 5)- 

r,5 
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where 
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r+S 



& = J- o J A ' dl > (2-4) 

and S = ±x, ±y (the lattice constant is taken to be unity for convenience). The 'on-site' and nearest-neighbour gap 
functions are defined as follows: 

A (r)=Vb< Ci (r)c T (r)); (2.5) 

A tf (r)=Vfr(c;(r)c t (r + ^) 

= -V s (c^r)c L (r + 6)), (2.6) 

assuming the existence oi- pairing in only the spin-singlet channel, in accordance with experimental results for the 
cuprate superconductors.E2l The mean-field Hamiltonian for the AvH model is written: 

r,5n,(T r,5 2 o,a" r ' a 

-^E ( r ) c t ( r + + H - c -] > ( 2J ) 

r,<? 

where r = mri + nr2, such that f\ = x + y and fi = x — y are primitive vectors of a single sublattice, and each lattice 
site has the two-point basis 0, x. Then, Sn = ±fi,±r2, <^20 = ±(^i + 7 ; 2),±(? : i — ^2), and 5 = ±x,±y. Throughout 
the remainder of this section calculations will be presented within the context of the EH model. Comparison of the 
above Hamiltonians indicates that analogous results for the AvH model can be obtained at any stage by eliminating 
the on-site gap function, reversing the sign of the kinetic term (keeping in mind the holes hop along second and 
third-neighbour links), and setting Vs = V/2. ■— ■ 
The Gor'kov equations can then be derived in the standard mannerEil 

e(r,r^c,o-e o (r,r^ WI o+E^ ( r I r ^^)[A (r'o-E A ^ r '')^( r '')]^ t ( r '^ r ^ w »); ( 2 - 8 ) 



J*(r, r', Wn ) = - G V - r, [A:(r") - £ W)M*")] S{i",i / , u n ), 



(2.9) 



where P$(r)X(r) = X(r + S) is the kinetic energy operator, Q and are finite-temperature single-particle and 
anomalous Green functions respectively, and Q° is the normal-state Green function in the presence of the external 
field. The Matsubara frequencies are oj n = itT(2ji + 1). The anomalous Green function T> is related to the gap 
functions: 



At(r)=TV J2^(r,r,u; n ), 
A* s (r)=TVs^2^(v,v + 5,u n ), 



(2.10) 
(2.11) 



Arbitrarily close to the s uper conducting critical temperature T c , the ratios A a /T c < 1 or identically zero (a — 0, 5). 
I terat ing Eqs. ( |2.§D and (2.9) up to third order in the gap functions, and making use of the conditions ( 2.1 0| ) and 
( [2.11 ), the self-consistent equations for the gap functions are immediately obtained: 

±A*(r) = - £ G°{r",T,-w n )A*{r")g°{r",r + a,u n ) 

r" ,u) n 

+ Yl 0V>>-w»)AVOSV>i>Wn)A(ri)S^ (2.12) 

r" ,ri ,r2 ,u3 n 
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where 



A*( X )=A:(x)-^AKx)P s (x), 



(2.13) 



and a = 0, ±x, or ±y. 

In the strong type-II limit, appropriate for the high-T c oxides, the penetration depth A(T) exceeds the coherence 
length £(T) and all other length scales. The single-particle Green function is then approximately translationally 
invariant: 



g°(v, r', Un ) « g°(v - r', u n )e~ 1 ^ £ A di 
where Q°(r — r',oj n ) is the normal-state lattice Green function in the absence of an external field: 

£°(x,W„) = V- — , 

^ lUJ n - £_k 



(2.14) 



(2.15) 



and the sum is over wavevectors in the first Brillouin zone. The dispersion relations for the EH and AvH models are 
respectively 



£™ = — 2t(cosk x + cosky) — [i 
£fc vH = 2tn(cosfci + cosfc2) + 4i 2 ocosfcicosfc 2 — fx, 



(2.16) 
(2.17) 



where k\ and fc 2 are reciprocal vectors of a given sublattice. Assuming the gap functions vary slowly compared 
with the characteristic length scale fop of the single-particle Green function (2.15), Eq. ( 2.1 2j ) can be expanded up 
to fourth-order in lattice derivatives. The justification for keeping higher-order gradient terms will be addressed in 
section V. The GL equations for the gap functions in the EH model can then be written: 

A* (r) = -TV a E Q ° ( mi + -^) g ° ( mi + ny + d'-a, W „)e"' i ^ A(rH5 " <3 ' ) 

W n a f 771,71 

(l - (in,) 2 - (yllyf - (xU x f - e™ (yH,) 4 - £ ™ 2 (xll x f (yH,) 2 W,(r) 



e ik-3 e -i|^A(r)-5 

K + 4 2 ) 2 



A D (r)|*A*(r) - 2|A (r)| 2 M(r) + 2|M(r)| :2 A*(r) 



-A* 2 (r)Af* (r) + A D (r)M 2 (r) - \M (r)| 2 M(r) \, 



where 



i {m z ) A_ z (r) = e l ^ A > {r) A z (r) - A_ z (r) 



47T 

zA z ) +i— A z (r) 



(2.18) 

(2.19) 
(2.20) 



with the forward difference operator defined by 



and 



zA z ) A_ 2 (r) = A_^(r + z) - A_ z (r) 
= A,(r) - A_ a (r), 



M(r) = ^ A|,(r)e ik - ? e% A(r) ' ? '. 

(5'— ±x,±y 



(2.21) 



(2.22) 



The coefficients of the gradient terms have been de fined in the Appendix for clarity. Note that all gap functions are 
taken at the same point r in the last term of Eq. ( [2.18| ); this condition will be relaxed in Section V. The gradient 
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terms with odd powers of xll x and/or yll y vanish for tetragonal systems. The GL equations for the gap functions in 
the AvH model are given by: 



A K r ) = ~FT / ] /] /] G° (mr i +nf 2 ,-uJn)G (mf 1 +nr 2 +5' -(5,w„)e~ 1 ^ 



*SfA(r)-(«-0 



u n c/ m.n 



1 - 4f (in,) 2 - e^f (yUy) 2 ~ 4f O^,) 4 - e^f {yli y f - (in,) 2 (yll y ) 2 AJ,(r) 



TV x ^ 



2 & (w - + ^ )2 



|M(r)| 2 M(r), 



where we have used 



(r x n ri ) = (in,) {yu y ) + (in,) + { y n y ) ■ 
(f 2 n r2 ) = - (in,) (yn y ) + (in,) - ( y n y ) , 



and the coefficients e are given in the Appendix. 
The current operator for the EH model is 

8, a 



-2tTVni(r')£(r,r>„)| , in^(r)/(r) = e i ^o A ^^f(r + 5) — f(r) 



2n 



n'(r) = -iV+ — A(r), 



(2.23) 



(2.24) 
(2.25) 



(2.26) 
(2.27) 



where ( 2.27 ) represents the continuum limit of the lattice curren t. Ite rating (2.8) and (2.9) to second order in gap 
functions and making the same approximations employed above, (2.26) becomes: 

j EH (r)=2tT Yl 



i l7r A < r K< s 2-<3i) 



e *o 



(-ljMt-ljWa^-Ri, Wn )a°(Ra - Ri - 3 X , -w„)A_ 5 (R 2 )S°(R 2 + a 2 , 



Ql ,Q2 

R,i ,R,2 



{ A; 2 (r) [m (in,)* + r> (yR,)*] A Ql (r) + A Q1 (r) [m! (in,) + n' (yll y )] A* 2 (r) } 



(2.28) 



where Ri = mx + ny and R 2 = mix + n'y. A similar expression for the AvH model can be obtained by setting 
t — > —tn, — 120 with (5 — > <5ii,<5 2 o and Ri = mfi + nf 2 , R 2 = m'f\ + n'r 2 , and making the replacements (in,) — > 
(in,) + (yU y ) and (yll y ) -» (in,) - (yRj. 

The various i ntegr als and sums appearing in the Ginzburg-Landau equations for the gaps functions (2.18), ( 2.23[ ) , 
and the current (2.28) can, in general, be determined only numerically. There is, however, one case which is analytically 
tractable: the EH model at low electron densities and weak to intermediate coupling. For V$ — and Vq — > — Vo, this 
limit would correspond to ordinary BCS theory. 



III. DETERMINATION OF T c 



At temperatures sufficiently near T c , we c an lin earize the gap equations ( 2.1 8| ) and ( 2.23 ). Making use of the 
definition of the normal-state Green function ( |2.15| ), we obtain for the EH model: 



A:(r) = TV, 



2 afc A;(r)-A;(r) 

k "i+e k 



TVi v a k [2o fc A;(r)-A;(r)] 

2 2.^ 



ASM = TVi E blA * dir) 



,,2 it!' 
k W n + ?fc 



(3.1) 

(3.2) 

(3.3) 
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where au = cosk x + cosk y and bk — cosk x — cosfcj,_-The s-wave and d-wave gap functions are related to the bond gap 
functions through the gauge-invariant definitions^ 



Mr) 
A d (r) 



;4(r), 



A x (r) + e l ^ (r) A_ x (r) + e-^ A » (r >A y (r) + {r) A_ y (r) 

?£ A * (r) A x (r) + e'^'WA^fr) - e^lS^W A„(r) - e^ (r) A_ y (r)' 



The equations that determine the s-wave and d-wave transition temperatures T s and 7^ are immediately 

\Vih{T s ) - 2] [V I (T S ) + 2] = VbVi7?(T a ); 



(3.4) 

(3.5) 
(3.6) 



(3.7) 



where 



h{T d ) 



h{T d ) = -, 
vi 



^ & V2T S 



, n = 0,l,2; 



tanh 



(3.8) 

(3.9) 
(3.10) 



It is clear from Eq. (3.7) that if V\ = 0, no s-wave superconducting instability can occur for positive temperature. 
The corresponding equations for the AvH model are 



V \ -» 1 + cos/^cosfcj + cosfci + cosfc 2 / £t 
1 = — ' tanh 



2T, 



V \ -> 1 + cos/cicosfc2 — cosfci — COS/C2 / £,k 

1 = tL : tanh 
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2T d 



(3.11) 
(3.12) 



While these equations are not analytically tractable, they have been shown numerically^ to strongly favour a d-wave 
transition temperature for all relevant hole concentrations. Optimal doping at approximately 25% filling (in the hole 
representation, where (lit) is calculated at T c ) yields a value of T d ~ 100 K. 

In order to analytically solve the sums 7i_ 4 , Eqs. (U) and ( |3~To| ), we make the standard transformation to an 
integration over energies, making use of the single-particle density of states (DOS) 



N(e) = J2S(e-e k ) 



JiV(0)+iV(±)hi|f |, ii\e\<D; 
\ 0, otherwise, 



(3.13) 



where D = At is the half-bandwidth and = + fJ-. The DOS is approximated by a constant plus a term reflecting 
the van Hove singularity at half-filling (e — 0), as is shown in Fig. |l]. The best fit is obtained when iV(0) = 0.31/1? 
(the DOS for ftee electrons in two dimensions is 1/ttD) and N{1) = 0.19/7* (note that the DOS at half-filling is 
approximately! 2 !] (2/7r 2 7>)ln|7>/£|). Using the relation afe = — 2(£ + /i)/7>, the sums 7q_2 can be solved to yield: 



I {T s )^2N'\n 



2D 



-N(l); 



(3.14) 



In 



- 1 



4m 
D 



D 



AT(1); 



(3.15) 



h{T s 



D 2 \T S 



4(7> 2 - 3^ 2 ) 
D 2 



N(0) 



2 

D 2 



ifiD + 6/i 2 ln 



D 



5/z 2 + 7> z 



N(l) 



(3.16) 
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where 



N' = N(0) - N(l)\n 



and 



For all terms proportional to N(\), we assume that T c <C D (weak-coupling), so that 

: ' ~ M \ _ f -1,£ < H; 



tanh 



2T 



{;!: 



(3.17) 



(3.18) 



(3.19) 



Note that in this lattice model, all interactions are instantaneous and therefore the bandwidth is the only possible 
energy cutoff. Since bk cannot be written exactly in terms of some simplifying assumption must be made in 
order to evaluate 1%. At low densities (the continuum limit of the lattice model), ~ tk 2 — D — fj, and thus 
b k « -2cos26>(£ + D + fi)/D. Then, 



D 2 



(D 2 - 4liD - 3^ 2 ) iV(0) 



2 

Z)2 



-2D 



(jjl + Df + ^(3^ + 4D)ln 



L» 2 - 8^L> - 5^ 2 



iV(l). 



(3.20) 



The equation for the s-wave transition temperature T s resulting from the application of Eq. (3.7) is, for Vq = 

-D 2 - Vi [6(i 2 N' - 2D 2 N(0) - (D 2 - 4fiD - 5/i 2 ) N(l)] \ 



T s = T*exp • 



(3.21) 



The corresponding equation for the d-wave transition temperature Td from Eq. (3.8) is 



T d = T*exp 



-D 2 - V 1 



H(4D + 3(i)N' - D 2 N(0) + ± (4D 3 + 7D 2 fi + UDfj? + 5^ 3 ) N{1) 



2((i + D) 2 ViN' 



(3.22) 



The transition temperatures for the cases V\ — t and V\ — it are shown as functions of chemical potential in Figs. |2| 
and U respectively. Near the bottom of the tight-binding band (fx « — D), an s-wave transition is strongly favoured for 
any Vq and V\, whereas a d-wave transition is favoured near half- filling (/i«0). The value of the chemical potential 
at which the preferred symmetry of the dominant gap function changes is extremely sensitive to the strengths of 
the respective interactions. As the on-site repulsion is increased, the magnitude of T s is suppressed while T4 is 
unaffected.Ej As a result, d-wave superconductivity is favoured for virtually all densities for sufficiently large Vq. It 
should be noted that the magnitude of the subcritical transition temperature, associated with the subdominant gap 
function, is in fact further decreased due to the presence of the dominant gap function. The relevant equation for the 



subcritical transition temperature, one of Eqs. (3.1)-(3.3), should have replaced by \J £ 2 



A° OM2 , where A° OM 

is the magnitude of the dominant gap funcion. Thus, the subcritical transition temperature, which is of dubious 
relevance in any case, can be taken to be identically zero without loss of generality. 



IV. CALCULATION OF THE GINZBURG-LANDAU COEFFICIENTS 



A. Extended Hubbard Model 



For sufficiently large systems, the lattice sums in Eq. ( 2.18| ) can be tra nsfo rmed into k-s pace integrals. Neglecting 
fourth-order derivatives and making use o f the expressions for the s-wave (3.4) and d-wave ( |3.5| ) gap functions as well 
as the normal-state Green function ( [2.15 ), the following three gap equations for the EH model are obtained: 
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d 2 k ( s* (r,k) 



d 2 z- 



d 2 k' 



i(k+k')-z 



(27r) 2 {-iu) n - £k'){iu n - £k) 

a: - 2TVt J (02 a fe{ as above } 



16|s(r,k)|V(r,k) +32|d(r,k)|V(r,k) + 16d* 2 (r, k)s(r, k) 
(cos 2 6»n 2 + sin 2 6»n 2 ) [s*(r, k) + d* (r, k)] 



(4.1) 
(4.2) 



A* d = 2TV 1 J2 



d 2 k fd*(r,k) 



(2^" K ^2 +C 2 ( w 2 +e 2 )2 



16|d(r, k)| 2 cT(r, k) + 32|s(r, k)| 2 cf(r, k) + 16s* 2 (r, k)d(r, k) 



dV 



d 2 fc' 



i(k+k')-z 



(27r) 2 (-icj„ - £k')(iu n - £fc) 



cos 2 6»n 2 + sin 2 6»n 2 ) [s*(r, k) + d*(r, k)] 



(4.3) 



where s(r, k) = ^A s (r) — jA (r), d(r, k) = ^■Ad(r), 9 is the angle between r and z (= mx + ny), and the continuum 
limit of ( gig ) 



n. 



+ — 4* W 



(4.4) 



acts only on the center of mass coordinate r. Eqs. ( 4.1 )-( 4.2 ) can be simplified by making use of the fact that 
A (r) = eA s (r), where e depends on temperature and chemical potential. This relation follows from the observation 
that ([O]) and (4.2) are not linearly independent, sinc e A (r) and A s (r) have the same symmetry. Near T c the 
magnit ude of the s-wave component is small; from we infer e rs VqI\/{1 + VqIq /2), w here Ip an d I\ are defined 
in Eq. ( |3.9[ ) and are evaluated in the weak-coupling and low-density limit in Eqs. ( |3.14 ) and (3.15). The integrals 
involving gradients of the gap functions can be greatly simplified by means of the identity: 



d 2 zx 2 G>,- W „)G°(z + o>„) = J J^L 
and similarly for integrals involving y 2 . Thus, 



v fex G°(k, w „; 



{2n) 2 oj 2 n +i 2 K 

^- I [l6a' fe 3 |A s (r)| 2 A:(r) + 846 2 |A (i (r)| 2 A:(r)+446 2 Af(r)A s (r) 



D 2 
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sin^n 2 + sin 2 fc,n 2 [2a' k A* s {v) + b k A* d (r)] 



(4.5) 



(4.6) 



^-5 k|A d (r)| 2 AS(r) + 8a' k 2 b k \A s (r)\ 2 A d (r) + Aa' 2 b k A* 2 {v)A d (v) 
(sin 2 ^!! 2 + sin 2 fc,n 2 ) [2a' fc A s * (r) + b k A* d (r) 



32 



where a' k = f - f . 

The appropriate free energy for a d-wave (or extended s-wave) superconductor on a tetragonal lattice is:tZI 
F s = F n + a s \A s (r)\ 2 + a rf |A d (r)| 2 + ^|A s (r)| 4 + /3 2 |A d (r)| 4 
+ /3 3 |A s (r)| 2 |A d (r)| 2 + 4 (A:(r) 2 A 2 (r) + A^(r) 2 A 2 (r)) 

+ 7s |IIA s (r)| 2 + 7d |IIA d (r)| 2 + 7 J(n j/ A s (r))*(n w A d (r)) - (U x A s (r)Y (U x A d {r)) + c.c. 



(4.7) 



(4.8) 
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where higher-order terms 



F 1 s = ?7s |A s (r)nA s (r)| 2 + 7s+ |n 2 A s (r)| 2 + 7s _|(n2-n|)A s (r)| 2 , 



(4.9) 



Ft = % |A d (r)nA d (r)| 2 + 7d+ |n 2 A d (r)| 2 + 7d _|(IT 2 - n 2 )A d (r)| 5 



(4.10) 



will be discussed in detail in Section V. In the present work, the s- wave and d-wave order parameters are respectively 
the gap functions A s (r) and A<j(r). The GL free energy ( |4.8| ) implies an s-wave (d-wave) component is induced 
whenever there exist spatial variations of the dominant d-wave (s-wave) order parameter. The magnitude of the 
'subdominant' order par ame ter is evid ently proportional to the coefficient 7 „ of the mixed gradient term. Minimizing 
and comparing to Eqs. (4J3) and (4/7) immediately yields: 



1 6 IT 



d 2 k 



d 2 k 



(2^ ui+e k 



(4.11) 



dPk {32< 4 , 2b\ , 32a' k 2 b 2 , 8a' k 2 b 2 } 

(2tt) 2 (cui+e k 



I) 2 



(4.12) 



{is , Ju , Id} = / (^2 Shl2 ^ 



{Aa' 2 , ~2a' k b k , b 2 } 



(u 2 n +e k ) 



2\2 



(4.13) 



The current den sity in the xy plane can be obtained either by evaluating (2.28) in the continuum limit or by minimizing 
the free energy (4.8) with respect to the vector potential: 



Aire 



Id 



AS(r)ffA d (r)+ 7s A:(r)rlA s (r)-i > (A:(r)n :c A d (r)+Aa(r)n a: A s (r)) 
+y>(A:(r)n y A rf (r) + A^(r)IL, A s (r)) + c.c. 



(4.14) 



where the factor 8 = (layer thickness/layer spacing) is introduced in order to model the layered structure of the 
cuprate superconductors. In th e pre se nt tw o-dimensional model, 8 — > 0, reflecting the inability of the system to 
sustain screening currents. Eqs. (4.11)-(4.14) are subject to the following boundary conditions: 



n • [ 7d nA d (r) + 7 „ {yU y - xU x ) A s (r)] = 0, 
n • [ 7s fiA s (r) + 7 „ (yU y - xU x ) A d (r)] = 0, 



(4.15) 
(4.16) 



where n is the unit vector normal to the surface of the superconductor. 

The coefficients of the GL free energy can be calculated analytically in the same low-density and weak-coupling 
limit employed earlier by setting a k « —2/i/D, b k w —2(fi + D)cos26 / D and sin 2 ^ ss 4(/x + D)cos 2 8/ D. Making use 
of 



d 2 k 



1 



7((3)N' 

(2n) 2 KT^F^'S^T^' 



(4.17) 



we immediately obtain: 



1 (l X 



2 U- L 

Vi \D A) 1 T s 



1_ f n + D 
Vi V D 



T 
T~d 



(4.18) 



A*, = *M (jj + 4 ) , - 



3 fn + D 



D 



f-+ £ V 

\D A) 



H + D 
D 



2\D AJ 



fi + D 
D 



, (4.19) 
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{7 S > Ju , Id} — 7^ j 2 (J^j + — 



D 



D 



D 



H + D 
D 



(4.20) 



where 



7C(3)A^ 

7T 2 T 2 ' 



(4.21) 



and vf. = D kp/4 « £>(/x + £>). From (|3.14J) and (J3.15Q we obtain: 



4 M 



D VoN'\n{^) + \ 



(4.22) 



The expressions for a s and ay are only valid near their respective critical temperatures. Since the subcritical tran- 
sition temperature is much lower than the dominant transition temperature, the corresponding coefficient for the 
subdominant order parameter can b e assu med approximately 1/V\ for all values of T < T c (i.e. the contribution of 
the appropriate integral in cither of ( 4.11 ) will be neg ligible). Note that in the limit V\ = 0, Vq — > — V q, A s and A d 
vanish; the relevant gap equation is therefore Eq. (4.1) with s(r, k) = Ao(r) and d(r, k) = 0. Then Eq. (4.8) becomes 
the two-dimensional continuum BCS GL free energy: 



F. 



BCS 



1 



T 

fi, 



7C(3)^V(0) 
16n 2 T 2 



Ao(r)| 4 



7C(3)^V(0) 
32tt 2 T 2 ' 



|nA (r)| 



/i 2 
8tt' 



(4.23) 



The GL coefficients (|4.18D - (|4.20D of the free energy (g 



imply that at the bottom of the band (/i = —D) the 
d-wave component vanishes, yielding a pure s-wave superconductor, while at half-filling (/i = 0) only the d-wave 
component remains (with the caveat that these analytical results become decreasingly valid near half- filling) . At 
densities intermediate between these two extremes all the GL coefficients are non-zero, leading to coexisting d-wave and 
s-wave order paramaters for any temperature T <T C and fiaite.. external magnetic field. For type-II superconductors 
described by ( |4.8| ) in fields just above H& , it has been foundEJ : E3E3 that the subdominant order parameter is nucleated 
in the vicinity of a magnetic vortex core whenever the coefficient of the mixed-gradient term 7^ is non-zero. Moving 
in the x = y direction from the center of the vortex, the induced subdominant order parameter reaches a maximum, 
then decreases algebraically; there exist extra nodes in the x = or y = directions. The maximum amplitude 
attained by the subdominant order parameter nucleated in the vortex core is given byiij 



A,(A rf ) max 
A d (A s )t^ 



\<*d(> 



16 ld(s 



a. 



(d) 



T 



T,, 



lis 



d{s) J ld(s) 



l-d(s) 



> T, 



s(d), 



(4.24) 



where A" lax is the maximum value of the induced subdominant s-wave order parameter compared wi th th e magnitude 
of the critical d-wave order parameter A^ ulk far from the vortex core, and vice versa. The estimates ( 4.24 ) are reliable 
as long as A s (A rf ) max < A,i(A s ) bulk /4, which is always the case sufficiently near T c : since A d (A s ) bulk 
the above equations indicate that 



A s (A d y 



{1-T/T C 



,3/2 



y/(l-T/T c ), 
(4.25) 



The subdominant order parameter must decay more rapidly near the transition temperature since it is induced through 
spatial v ariati ons of the critical order parameter. 
From ( 4.20 ), we have: 



Id 

Is 



M + De/4 



H + D 



fi + De/4 



T s > T d , 



(4.26) 
(4.27) 



The gradient coefficient ratios ( [l.2q ) and (4.27), which govern the magnitude of the subdominant order, are compared 
with the appropriate numerical results in Fig. [4] for the special case Vq — 0. The analytical results capture the essential 
physics in their regime of applicability, i.e. at low densities and weak to intermediate coupling (better quantitative 
agreement at intermediate coupling can be obtained by including in ajfc, bk, etc. terms higher order in £&). Near 
half-filling (zero- filling), where d-wave (s-wave) superconductivity is most stable, 7^/7^ {^ v /ls) grows with increased 
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coupling. As a result, a significant component of the subdominant order parameter would only exist for stronger 
coupling in these two density regimes. It should be kept in mind, however, that as the on-site repulsion Vq (and 
therefore e) increases from zero, the magnitude of the s-wave component is suppressed for all densities; a d-wave 
subcomponent would be thereby enhanced in a bulk s-wave superconductor. The magnitude of the subdominant 
order parameter can be significantly larger even for weak coupling at intermediate densities, however; not only is the 
gradient coefficient ratio enhanced, but also the ratio | ctdom | /c*sub since a su b as the subdominant superconducting 
instability is approached. (Fig. |^ in the following section provides further quantitative details). 

There are in general two characteristic length scales £ s and ^ governing spatial variations of the s-wave and d-wave 
order parameters in the vicinity of a vortex core, respectively. Near T Cl however, the induced subdominant order 
parameter is negligible compared to the critical order parameter by ( 4.25| ). Keeping only terms in the free energy 



(4.8) involving the dominant order parameter, the superconducting coherence lengths £(T) and penetration depths 



A(T) near T c can be crudely estimated: 



Ud)(T) = (4.28) 



A s(rf) (W .. 4 y i( f (4.29) 

V <H 47r ) 7s(d)K(d)l 

where the flux quantum 4>q = he/ e and both lengths are given in units of the lattice constant a. Inserting the relevant 
nalytical expressions from Eqs. Q4.18 )-(4.20), valid in the low-density, weak-intermediate coupling limit, we obtain: 



a. 



UT) ~\2^f)p(l-T/T s ) ' U{T} ~ \fif) V 4(1 -T/T d )' (4 ' 30) 



A S (T) 



n 3 5v 2 Jl 



T/T s ) 



4n 3 Sv 2 F (l 



T/T d 



(4.31) 



Note that the penetration depths are effectively infinite for an isolated layer. The magnitudes of the coherence lengths 
can be estimated for parameters stabilizing either s-wave or d-wave superconductivity (depending on the magnitude 
of Vq). For example, with fj, = —2.0t, T c ~ 0.2t (Vi ~ 3i), and T = 0.9T C , we obtain ^ ~ 6a and £ s ~ 8a for d-wave 
and s-wave superconductors respectively. The coherence lengths tend to become progressively shorter with increased 
electron density (Fig. |^ in the next section gi ves f urther details) . 



It is useful to compare the GL free energy (4.8) with that derived by Ren et 



Defining s = — (-^ 



A., and 



d — Ad (we drop reference to the center of mass coordinate r for convenience), we obtain: 



F = F 



T 

fi, 



T 
fi. 



|d| 2 +4 7 {| 



r |d| 4 + 2| S | 2 |d| 2 



)} 



+ 7 



(U y s)*(llyd) - (U x s)*(n x d) +c.c.]} 



8tt' 



(4.32) 



While the various coefficients appearing in the above free energy appear t o be similar to those found by Ren et 
at, one important distinction must be emphasized. The above free energy (4.32) mixes a d-wave order parameter 
with an extended s-wave order parameter. The analytical results suggest that only near half-filling (fi — > 0) does 
the contribution of A s to s vanish. At these densities, however, the approximations employed in the microscopic 
calculation above are not strictly valid; indeed, the numerical results clearly indicate the existence of a non- vanishing 
s-wave component for T < Td even in the absence of an on-site interaction (see Fig. ^) . For any other density (except 
fi — > —D, where the d-wave component vanishes) at temperatures below T c , both A s and will be non-zero in 
the vortex core since both are generated by the same nearest- neighbour pairing interaction V\ . There is therefore no 
limit in which the present model reduces to that of Ren et at, i.e. where an isotropic s-wave order parameter derived 
from a repulsive on-site interaction is alone nucleated by spatial variations of a dominant d-wave order parameter. 
The technical difficulty which necessitated the implementation of the Pade approximation in Ref. |l^ was due to the 
indistinguishability of A and A s in the continuum. This problem can always be avoided by starting with a lattice 
model which gives rise to a d-wave order parameter, then taking the appropriate continuum limit. 
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B. Antiferromagnetic-van Hove Model 



Neglecting fou rth-o rder derivatives, the equations for the gap functions in the thermodynamic limit of the AvH 
model from Eq. ( 2.23| ) are: 



TV 



E 



d l z- 



IX. v 



d 2 z — 



d 2 k 1 


(c k A* s (r) 


1 




(2tt) 2 1 




4K + 




' d 2 k' 




-k')-z 




(27r) 2 (-iu> n - £ k 


)(iw„ - 




d 2 k i 




1 




(2tt) 2 1 






a 2 


' d 2 k' 


e *(kH 


-k')-z 





c 2 |A s (r)| 2 A;(r) + 2c fc d fc |A d (r)| 2 A:(r) + c k d k A* 2 (r)A a (r) 



n 2 + sin26» (n 2 - n 2 ) 



[c fe A:(r)+e fc A;(r)] ; 



(4.33) 



(27r) 2 (-zw„ - £ k >){i<JJ n - 6c) 



d 2 |A d (r)| 2 AS(r) + 2c fc d fe |A 5 .(r)| 2 A2(r) + c k d k A* 2 {v)A d {v) 
n 2 + sin20 (n 2 - H 2 )] [e k A* s (r) + d k A* d (r)} ), 



such that 

c k = 1 + cosfciCOsfc2 + cos/ci + cosfc2; 
rffe = 1 + cosfcicosfc2 — cos/ci — COS/C2; 
efc = sinfcisinfc2. 

From these, and making use of (|4.5|), we can immediately obtain the appropriate free energy (fbq) where 



r I I [ d 2 k {C k , d fc } /a 

la s , ay^ = / - — ttt tanh — 

1 ' S V A J (2^) 2 6 V2T 



(4.34) 



(4.35) 
(4.36) 
(4.37) 



(4.38) 



{01 ,02,03, 04} 



T 



d 2 k {c\ , d\ , Ac k d k , c k d k } 



16 As J {2n) 2 



2\2 



(4.39) 



{7s , Id} = T 



d 2 k sin 2 /ci {t\\ + 2i2o c os/c2) 2 + sin 2 /c2 (in + 2^2ocosfci) 2 



{cfc , d k }, 



(4.40) 



7^ 



-2T 



As 



d 2 k e 2 (in + 2t 2 ocos/ci) (in + 2i 2 ocos& 2 ) 
(2^ 



(4.41) 



Note that /33 = 4/34, as was found for the EH model ( 4.32 ) and in Rcf. |l9j. This generic feature of the free energy (4.S) 
is a direct consequence of the symmetry between the s-wave and d-wave order parameters and the underlying bond 
(or directionally-dependent) gap functions from which they are derived. Due to the complicated angular dependence 
of the AvH dispersion ( [2.17 ), the above expressions above are not analytically tractable, howe ver. At optimal hole 
doping ((rihoio) = 1 — ("-electron) ~ 0.2), the d-wave transition temperature found using Eq. (3.12) is ~ 100 K, and the 
coefficients of the GL free energy have been evaluated numerically: 



T 

T, 



F s = F n - 11.3 I 1 - — I |A d | 2 + 10.8IAJ 2 + 38.7IAJ 4 + 5720|A d | 4 + 603 



A.,I 2 IA 



2.69IIIAJ 2 + 7.38|nA (i | 2 + 1.50 



(n^nn^) - (n x A s )*(n,A d ) + c.c 



li 

87T' 



A! 2 A 2 + Af A 2 



(4.42) 



where the coefficients 0, 7 have been evaluated at Tj. 

The maximum value of the s-wave component nucleated in the vortex core, calculated numerically using (4.24), 
(4.11), (4.13), (4.38), and (4.40), is shown as a function of hole density in Fig. ||. Evidently, in the AvH model the 
induced s-wave component can be at most 5% of the bulk d-wave value (corresponding to (n/j) ~ 25% and T — ► 0, 
where it is assumed that aa{T) oc (1 — T/TA ; VT). This is in spite of a large temperature-independent gradient 
coefficient ratio (7^/7^ ~ 20% near optimal doping) that is found to increase monotonically with doping up to the 
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largest densities studied. By contrast, with a suitable choice of parameters (low T,Vo, large (n^), Vx), the EH model 
can allow for s ma x/dbuik ~ 20 — 30%. As the on-site repulsion is increased, however, the induced s-wave component 
is suppressed. The reduction in the magnitude of the s-wave component in the 'overdoped' AvH model may be 
partly due to effective on-site repulsion, built into the Hamiltonian by constraining holes to move within a single spin 
sublattice of an antiferromagnetic background. 

The GL d-wave coherence length £(T), calculated numerically with the aid of ( 4.28 ), is shown as a function of 
hole concentration for T = 0.9Td in Fig. |(| As expected, £(T) becomes shorter with decreasing (rih). For hole 
densities up to 30% filling, the-AvH model leads to £ ~ a/ y/1 — T/T d where a is a lattice spacing, in agreement with 
experiments for the cuprates.ES The EH model yields similar d-wave coherence lengths for V\ = 2t, 3t. These short 
£d strictly violate the initial GL condition that the gap functions be more slowly varying than the Fermi wavelength 
kp 1 ~ la — 2a, the characteri stic length scal y f|f ft hp fcrmionic excitations. The close agreement, however, between 



the results of integrating Eqs.pd4.8|) and ( 4T4 )tjE3E9 and numerical investigations of short-^ superconductors within 
Bogoliubov-De Gennes theorjJl3 indicates that the GL equations derived thus far capture at least some of the essential 
physics. 



V. EXTENSION OF THE GL EQUATIONS 



The current and magnetic field distribution around an isolated vortex of a d-wave superconductpij.is found to be 
fourfold symmetric whenever an s-wave component is induced in the vicinity of the vortex core.E3E3 It has been 
recently shown, how ever , that a similar anisotropy results from the inclusion of higher-order d-wave gradient terms in 
the GL free energy (L8) even in the absence of an s-wave component £3 The connection between these observations 
can be made apparent by integrating out of the free energy the degree of fr eedom associated with the subdominant 
order parameter. Given a bulk d-wave superconductor, for example, Eq. (4.8) can be minimized with respect to A*(r), 
yielding to lowest order in (1 — T/Td) 



A s (r) = -^(n 2 -H 2 )A d (r). 



(5.1) 



Note that s ince A d (r) varies on the length scal e £d , we immediately obtain A s ~ 7iVc* s £j ~ (1—T/Tj) 3 ' 2 in agreement 
with (4.25). Upon substitution of (5.1) into ( |4.8| ), the leading correction due to the induced s-wave component will 

2 

be of the form — ^|(n^ — H|)Ad| 2 . The free energy for an inhomogeneous d-wave superconductor then becomes: 



F s = F n + a d \A d (r)\ 2 + /3 2 |A d (r)| 4 + 7d |IIA d (r)| 2 
+ r /d |A rf (r)nA d (r)| 2 + 7d+ |II 2 A d (r)| 2 + 



Id- 

oc. 



(IE 



n 2 )A rf (r)| 2 + -, 



(5.2) 



which clearly indicates that the induced s-wave order parameter breaks rotational symmetry in precisely the same 
way as a fourth-order d-wave gradient term. 

In order to assess the relative importance of the induced subcritical component in eliciting anisotropy, the coefficients 
of the higher-order terms included in (4.8), and 7d( s )- in particular, must be determined microscopically The 
derivation of r\ s and rj d appearing in (4.9) and (4.10) requires extending the second term on the right hand side of Eq. 
(2.12) to include gradients of the gap functions. Applying similar techniques to those described in previous sections, 



one obtains for the EH model: 



(2ir) 2 r,.,2 _i_ L 




32tt 4 T 4 



D 



D 



(5.3) 



where the analytical solution is valid in the continuum limit, i.e. at low electron densities. For the AvH model, 



d 2 k €l 



(v fel ^) 2 + (v fe2 &) 2 



{4 , 4}. 



(5.4) 



The c oeffici ents of the fourth-order gradient terms can be derived by evaluating in the thermodynamic limit terms in 
(2.18) and (2.23) hitherto ignored. After some manipulation, one obtains 
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{ls± , ld±\ = > / d z / - — — — —a; (as ± 3y ) \ 4a k , b%\ , 

{%± , 7, ± } AvH = -jE/^/ • • ■ + ^) 2 [(-i + -2) 2 ± 3 (n - r 2 ) 2 ] {c fc , d k } , (5.5) 



where the ellipsis represents the /c-dependent part of the integrand. The fourth-order gradient terms can be evaluated 
analytically in the EH model for weak to intermediate coupling and low densities: 



b s+ , J d+ , 7- , 7,-} « ^ . £ + 7 = , ^ .0 ,~ ^ . (5.6) 



_ 31C(5)iV' f 3 / m £ 

1 ~ 256 VttT/ j 2 VD 4 




The analytical results obtained for the EH model indicate that all the higher-order terms in the free energy (4.G) 



and (4.1C) have a negative sign compared with the second-order terms. The overall sign of these higher-order terms 
is of no consequence, however, provided the order parameters are sufficiently small and slowly varying; while both 
conditions will be satisfied near T c , we also find r ) s (d)+lls(d) ~ (vf/kT) 2 — > as [i — > — D. At low densities, the EH 
model predicts a vanishing 7 S _ coefficient, in accordance with the expectation that the free energy of a continuum 
s-wave superconductor should have spherical symmetry. Furthermore, ^d+/ld- = 5/2, as was found by Ichioka et 
alB 

The analytical results for the EH model cannot adequately determine whether anisotropy in the free energy is 
primarily due to the existence of the subdominant order parameter or fourth-order gradients, however. For a bulk 
s-wave superconductor, the induced d-wave order parameter clearly breaks circular symmetry since 7„ ^ for /i > — D 
while 7 S _ = 0. Yet stronger coupling or lat tice effects could lead to a non-vanishing 7 S _ which could compete with 



the 7y/ci!(i coefficient in the free energy ( [5.2| ). Assuming that a s w 1/Vi, we find for a bulk d-wave superconductor 

-il 196C(3) 2 ViA^(0) / /i 



a sld . 31C(5) \D 



(5.7) 



which is of order —V\/t. This result indicates that the 'subcritical' coefficient 7^/a s and the 'asymmetric' coefficient 
7d_ give similar contributions to the anisotropy. It should be emphasized, however, that d-wave superconductivity is 
not favoured at densities for which these analytical results are strictly valid. 

The ratios ^lt/ad{ s )ls{d)~ have been calculated numerically for T = T c and are shown in Fig. |7]; virtually indistin- 
guishable results have been obtained for all temperatures 0.5T C < T <T C (not shown). The numerics make clear that 
the contribution to anisotropy of the subcritical coefficient is at most comparable to the asymmetric coefficient for 
most densities and coupling strengths. Furthermore, sinc e both coefficients have the same sign, their contributions 



to the asymmetric gradient term in the free energy (5^) are in fact competing. Only for densities very near the 
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crossover from bulk s-wave to d-wave condensation (or vice versa) in the EH model, or for large hole densities in the 
AvH model, does the subcritical coefficient dominate; in this regime a — > 0. 

It is not presently clear if the overall sign of the anisotropic fourth-order gradient term in the free energy ( |5.2| ) has 
any physical significance. Previous microscopic investigations of the EH model within the context of Bogoliubov-De 
Gennes theorytS demonstrated marked anisotropy in the structure of the critical d-wave component near the core 
of an isolated vortex. Parameters chosen correspond to a substantial s-wave component nucleated near the vorti 
core, i.e. 7 2 /a s 3> 7d-, and therefore a large negative coefficient for the asymmetric gradient term. Recent work,! 
however, indicates that this anisotropy persists even for densities approaching half-filling, where 7d_ dominates and 
the overall sign of the gradient term is positive. 

For completeness, the coefficients rj s (d) an d J s (d)+ have also been determined numerically for both the EH and AvH 
models. Both coefficients are always negative and in general it is found that for the EH model ?7s(d) ~ 107s(d)+ ~ 
100|7 s (d)— |, and 7 s (d)+ is of the same order as j s (d) at T c . For the AvH model at optimal doping and T c , we obtain: 

F s AvIi =F n - 11.3 (l - £-) |A rf (r)| 2 + 5720|A d (r)| 4 + 7.38|flA d (r)| 2 



- 16300| A d (r)nA d (r) | 2 - 16.5|II 2 A d (r)| 2 + 2.65|(n 2 - n 2 )A d (r)| 2 + (5.8) 

The relatively large values of T} s (a\ and "f s (d)+ compared with the magnitude of 7 s (d) clearly demonstrates that the 
GL theory derived herein is only strictly valid quite near T c . 
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VI. SUMMARY AND DISCUSSION 



The primary objective of the present work has been to derive the Ginzburg-Landau equations for a d-wave super- 
conductor using two microscopic lattice models which have been previously used to describe the high-T c oxides: the 
extended Hubbard model and the Antiferromagnetic-van Hove model. In so doing, it has been possible to quantita- 
tively investigate how the lattice and external magnetic field generate, and govern the interplay between, co-existing 
s-wave and d-wave order parameters. In addition, the relative magnitudes of the various GL coefficients, as well as 
their temperature and density dependence, have been ascertained. 

Phenomenological GL theory has enabled much progress to be made recently toward understanding the structures 
of isolated vortices and tlie.iiartex lattice for d-wave superconductors. Of particular current interest is the theoretical 
prediction by Franz et alE3t3 of an oblique structure for the vortex lattice near H C 2, i.e. the Abrikosov lattice would 
be intermediate between the usual triangle and a square. The degree of 'obliqueness' is mostly dependent on the 
gradient coefficient ratio 7^/7^. The coefficient 7„ governs the extent an s-wave component is induced by spatial 
variations of the dominant d-wave order parameter, and characterizes the degree of fourfold symmetry in the free 
energy. For 'Jv/jd = 0, the s-wave component vanishes, yielding a triangular lattice. The Abrikosov lattice deforms 
continuously away from a triangle as 7„/7d is increased; for 7^/7^ = 0.45, the angle between primitive vectors (j> = 76°. 
For "fu/ld ~ 0.6 and higher, the flux lattice is square. 

The present work demonstrates that microscopic models used to describe the high-T c oxides can predict a significant 
admixture of an s-wave order parameter in the mixed state of a d-wave superconductor. One consequence is the 
deviation of the flux lattice from that of a triangle. It has been found that within a broad and experimentally relevant 
parameter space both microscopic mndpls yield a gradient coefficient ratio 7^/7^ ~ 0.1 — 0.4. This is consistent with 
two recent experimental observationsE3c3 for YBCO of flux lattices with rj 73° and 77°. It is not yet clear, however, 
whether the a — b anisotropy, associated with the orthorhombicity of YBCO is alone sufficient to account for the 
distortion in the flux lattice.tZl 

It is presently uncertain whether the s-wave component nucleated in the vicinity of vortex cores is in fact required 
to induce a significant deviation from a triangular vortex lattice. It has been shown that a fourth-order gradient 
term in the d-wave order parameter also introduces fourfold symmetry into the free energy; over much of the phase 
diagram, the contribution of this term to anisotropy is comparable to that of the s-wave component. This important 
issue will addressed in future work. 

Wherever possible, comparison has been made with previous derivations of the GL coefficients within continuum 
models. It should be emphasized, however, that lattice models not only provide considerably more information 
regarding the density and coupling-dependence of s-wave and d-wave admixture in the vortex core, but also avoid the 
technical difficulties (i.e. the application of the Pade approximation) encountered in continuum models. In particular, 
the lattice models clearly indicate that both s-wave and d-wave components must always coexist in the vortex core 
for all temperatures below T c , regardless of the symmetry of the bulk order parameter. 
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APPENDIX: COEFFICIENTS OF THE GL GRADIENT TERMS 



The coefficients of the gradient terms appearing in the GL equations for the gap functions ( p.l§| ) and (2.23) are, 
for the EH model: 
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For the AvH model, we obtain: 
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FIG. 2. The transition temperatures T c are given as functions of chemical potential jj, for nearest-neighbour attraction V\=t 
and Vo = 0. The results for T s and Td are given by solid and dashed lines, respectively. The analytical results (lighter lines) 
are found to compare well with the numerical results (darker lines) in the applicable low-density regime. 
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FIG. 3. The transition temperatures T c are given as functions of chemical potential /i for Vi = 3i. The solid lines correspond 
to the numerical evaluation of T a for on-site repulsion Vo = 0, 2.5t, 3.5i, and ~ oo; linewidth decreases with increased Vb. 
The d-wave transition temperature Td (dashed lines) is unaffected by changes in Vb. Note that the analytical (lighter dashed 
line) and numerical (darker dashed line) results for T d still agree closely at low densities in this intermediate-coupling regime. 
The correspondence between the numerical and analytical results for T s (not shown) continues for intermediate coupling but 
is found to improve with increased Vb, or decreased T s . 
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FIG. 4. The ratio of the mixed gradient coefficient 7„ to the ordinary gradient coefficient of the dominant order parameter 
7s or "fd is given as a function of chemical potential /i. The solid and dashed lines correspond to numerical and analytical 
results respectively. The solid lines become progressively darker as Vi is decreased from 3t to t in increments of t; Vo is taken 
to be zero. The discontinuity reflects the transition from an s-wave to d-wave bulk superconductor; the dashed line assumes 
that this occurs for /i = — 1.65t. All values are determined at the appropriate T c . 
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FIG. 5. The relative magnitude of the maximum s-wave component induced in the vortex core s max /dbuik is calculated 
numerically and given as a function of hole density {nu) for temperatures T = 0.97^. The dashed line, corresponding to results 
for the AvH model, indicates that the induced s-wave component is largest near optimal doping ((rih) ~ 20%). Results for the 
EH model are given for comparison, where in the hole representation {rih} = corresponds to [i = 0; solid lines (in order of 
decreasing boldness) correspond to Vi = t, 2t, 3i with Vb = 0, while the dotted line is for V\ = 3i, Vo = 4t. 
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FIG. 6. The numerical d-wave Ginzburg-Landau coherence length £d(T), measured in units of a lattice spacing, is shown for 
T = 0.9Td as a function of hole concentration (n/t). The dashed line corresponds to the AvH model while the solid lines (in 
order of decreasing boldness) correspond to the EH model for Vi = t, 2t, 3i and Vb = 0. 
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FIG. 7. The ratio 7^/Qd(s)7s(d)- is shown as a function of electron concentration (n e ) for T — T c . Solid lines (in order of 
decreasing boldness) correspond to the EH model for Vi = 1.3t, 2t, 3t with VI) = 0, while the dotted line is for the EH model 
with Vi = 3t, Vo = 4t. The dashed line gives the results for the AvH model in the electron notation such that (n e ) = 1 — 
To the left (right) of the arrow is shown 7^/Qd7s- (j„ / a s Jd-) corresponding to bulk s-wave (d-wave) superconductivity. 
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